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Generally,  the  Permanent  Magnet  Tubular  Generators  (PMTG)  are  employed 
in  the  free-piston  generator  systems.  In  this  paper,  the  PMTG  with  the  eore 
made  of  Soft  Magnetic  Composite  (SMC)  is  employed,  and  then,  the  PMTG 
with  the  solid  SMC  core  has  cylindrical  symmetry  and  also  the  outer 
diameter  of  the  machine  reduces  in  comparison  to  the  ones  with  multipart 
laminated  steel  cores  in  which  each  part  is  radially  placed  along  the  air-gap 
and  there  is  unexploited  space  between  them.  Here,  a  Finite  Element  Method 
(FEM)  is  developed  in  which  the  reciprocating  motion  of  the  translator  is 
considered  by  moving  mesh  techniques  and  moreover,  the  Jiles-Atherton 
(JA)  scalar  hysteresis  model  is  incorporated  in  the  EE  model  to  take  into 
account  the  hysteresis  losses.  Then,  the  copper  and  core  losses  are  computed 
and  the  equivalent  circuit  of  the  generator  is  derived.  The  calculated  losses 
are  compared  with  the  ones  obtained  from  the  analytical  methods. 
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1.  INTRODUCTION 

In  recent  years,  the  importance  and  applications  of  linear  electric  machines  has  been  increased, 
drastically.  Similar  to  the  rotating  electric  machines,  the  linear  electric  machines  are  also  divided  into 
different  types  in  terms  of  principles  of  their  performance.  The  reciprocating  permanent  magnet  (PM)  linear 
electric  machines  have  been  found  the  more  importance  due  to  having  the  high  force  density,  high  efficiency, 
and  simple  control  configuration  [1],[2].  The  reciprocating  PMTG  in  which  the  crank  shaft  is  omitted,  is 
known  as  the  free-piston  generator.  For  this  reason,  the  mechanical  loss  is  low  in  the  reciprocating  PMTG. 
These  electric  machines  can  be  used  in  hybrid  electric  vehicles  [3]  and  for  electricity  generation  from  sea 
waves  [4], [5]. 

In  previous  researches,  the  more  attentions  have  been  paid  on  the  configurations  and  mechanical 
problems  of  free-piston  generators,  and  the  loss  calculation  has  not  received  much  attention  [6].  The  exact 
analytical  modeling  of  linear  electric  machines  is  laborious  due  to  their  complex  configuration  [7], [8].  For 
this  reason,  the  finite  element  method  (FEM)  is  used  in  this  paper  for  modeling  the  reciprocating  PMTG  with 
SMC  core.  The  aim  of  this  research  is  the  consideration  of  iron  loss  in  free-piston  generators  with  the  SMC 
core.  The  SMC  material  has  high  hysteresis  losses,  and  so  considering  the  hysteresis  loop  is  very  important 
for  modeling  the  free-piston  generators  with  the  SMC  core.  The  J-A  method  is  used  for  modeling  the 
hysteresis  loop.  In  this  research,  the  SMC  core  is  used  for  the  first  time,  and  with  respect  to  the  specific 
features  of  SMC  material,  in  addition  to  the  loss  calculation,  the  synchronous  inductance  is  calculated,  the 
equivalent  circuit  is  extracted,  and  finally  the  characteristic  of  voltage,  output  power  and  the  input  power  in 
terms  of  displacement  is  also  drawn. 
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The  COMSOL  software  package  is  used  for  simulation  by  FEM.  For  verification  the  parameters 
used  in  simulation,  the  analytically  calculated  loss  is  compared  with  the  corresponding  results  obtained 
through  FEM. 

This  paper  is  organized  as  follows:  the  configuration  of  reciprocating  PMTG,  and  the  parameters 
used  in  FEM  are  introduced  in  section  II.  The  techniques  for  loss  calculation  are  also  reviewed  in  section  II. 
The  results  obtained  through  FEM  are  presented  in  section  III.  The  losses  are  analytically  calculated  in 
section  IV.  The  synchronous  inductance  is  calculated  in  section  V  and  the  equivalent  circuit  is  then  drawn. 
Section  VI  gives  the  conclusions. 


2.  MACHINE  CONFIGURATION  AND  ANALYSIS  METHOD 

The  reciprocated  PMTG  analyzed  in  this  paper  was  introduced  in  [3].  This  generator  is  used  in  the 
combustion  motor  of  hybrid  vehicles  as  the  free-piston  generator.  Figure  1  shows  the  configuration  of  this 
generator.  The  winding  of  this  generator  includes  6  coils  connected  in  series,  and  its  translator  has  4  PM 
poles  with  alternative  polarities.  The  only  possibility  for  manufacturing  the  tubular  stator  core  with  laminated 
ferromagnetic  materials  is  shown  in  Figure  2  [9], [10].  In  this  research,  material  SMC-Somaloy500  is  used  to 
manufacture  the  solid  tubular  core.  SMC  has  specific  features  such  as  negligible  eddy-current  losses, 
flexibility  in  getting  form,  the  suitable  cost,  and  high  Curie  temperature  [11]. 


Figure  1 .  Reciprocating  PMTG  in  free-piston  generator  system 


Figure  2.  Laminated  axial  view  of  tubular  machines 


The  hysteresis  losses  of  translator  have  been  ignored.  The  specifications  of  this  machine  are 
presented  in  Table  1.  As  shown  in  Figure  3,  the  displacement  profile  of  translator  along  its  axis  is  sinusoidal 
with  angular  speed  150  rad/sec.  Figure  4  shows  the  B-H  curve  of  core  material  [12].  The  resistive  load  150  Q 
is  connected  to  terminal  of  generator. 
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Table  1.  The  PMTG  specifications 


Parameter 

Value 

Output  power 

116.9235  W 

Output  voltage 

102.0676  V 

Pole  number 

4 

Diameter 

135.5802  m 

Turn  number 

296 

PM  thickness 

18.415  mm 

PM  length 

37.719  mm 

Air-gap  between  adjacent  PMs 

4.191  mm 

Tooth  width 

20  mm 

Tooth  depth 

20.955  mm 

Slot  width 

20  mm 

Slot  depth 

20.955  mm 

Yoke  width 

14.7066  mm 

Core  length 

25 1 .46  mm 

Air-gap  length 

1.651  mm 

Shaft  diameter 

26.035  mm 

PM  remanence 

0.8  T 

Global: 


Time  (s) 


Figure  3.  Displacement  profile  of  translator 


BH(t) 


Figure  4.  Anhysteretic  Magnetization  curve  of  SMC  SomaloySOO 


2.1.  Modelling  method  and  hysteresis  loss  calculation 

The  most  accurate  method  for  hysteresis  loss  calculation  is  using  the  hysteresis  loop.  There  are  two 
existent  methods  of  preisach  [13]  and  Jiles-Atherton  (J-A)  [14]  for  considering  the  hysteresis  loop.  The  J-A 
model  is  more  feasible  and  simple  than  preisach  model  while  considering  the  sinusoidal  excitation  and  the 
SMC  material.  Under  this  condition,  the  computation  time  of  FEM  by  J-A  is  the  half  of  it  by  preisach  model 
[13].  Therefore,  the  J-A  method  is  used  in  this  paper  for  modeling  the  hysteresis  [15].  The  vector  equations 
of  J-A  method  are  as  follows: 


dt  dt  dt 


coth- 


He\  a 


a,  I I 


\H\ 


(1) 

(2) 
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H^=H+aM 


(3) 


dM,„ 

dt 


dj^ 

dt 


(4) 


M,„=c.(c-r.(M-MJ 


(5) 


M  ,  and  M^^^are  the  total  magnetization,  the  anhysteretic  magnetization,  the  irreversible 

magnetization  and  the  reversible  magnetization,  respectively.  Also,  H  and  represent  the  magnetic  field 
vector  and  the  effective  magnetic  field,  respectively.  Here,  COMSOL  Multiphysics  Software  is  employed  in 
which  the  corresponding  external  electrical  circuit  is  coupled  to  the  electromagnetic  equations  that  are 
formulated  based  on  the  magnetic  vector  potential,  and  regarding  the  implementation  of  the  above  mentioned 
equations,  more  details  can  be  found  in  [16].  Empirical  parameters  M^,a,  a,  k  and  c  are  the  parameters  of  the 
J-A  hysteresis  model,  whose  values  are  extracted  from  experiments  [17]. 


3.  FEM  RESULTS 

Figure  5  shows  the  refined  mesh  of  the  geometry.  The  flux  density  distribution  of  reciprocating 
PMTG  is  shown  in  Figure  6.  As  shown,  the  maximum  flux  density  occurs  when  the  translator  is  in  the  middle 
position,  and  the  minimum  flux  density  also  occurs  when  the  translator  is  in  the  end  positions.  Figure  7 
shows  the  flux  tube  distribution  in  the  core. 


Figure  5.  FEM  mesh  of  the  PMTG  mode 
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Figure  6.  Flux  density  distribution 
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Figure  7.  Magnetic  flux  tubes 
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The  hysteresis  loop  for  each  tooth  is  different  from  other  teeth.  Therefore,  the  hysteresis  loss  should 
be  calculated  for  each  tooth,  separately.  Since  the  magnetic  flux  density  in  yoke  is  affected  by  teeth,  so  the 
hysteresis  loop  should  be  extracted  for  each  part  of  yoke,  separately.  Figure  8  shows  the  hysteresis  loop  for 
middle  tooth.  The  hysteresis  loss  is  calculated  as  follows: 


xf 


(6) 


where  is  the  hysteresis  loss,  B  is  the  magnetic  flux  density,  is  the  volume  of  the  core,  and  f  =11 T 
is  the  frequency.  The  sum  of  area  of  the  hysteresis  loop  in  radial  and  vertical  directions  is  equal  to  the 
hysteresis  loss  energy.  For  the  analyzed  PMTG  in  this  paper  =  5.\ll{watf) . 
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Figure  8.  The  hysteresis  loop  in  the  middle  tooth  in  radial  direction 


The  core  of  stator  and  translator  should  be  manufactured  in  solid  form  due  to  the  configuration  of 
reciprocating  PMTG,  and  consequently  the  eddy-current  loss  is  high.  For  this  PMTG,  total  eddy  current  loss 
is  low  (0.562  watt)  due  to  using  the  SMC  material  in  stator  and  translator  cores.  The  eddy-current  loss  is 
calculated  as  follows: 

(7) 

a 

where  cr  is  the  conductivity  coefficient  of  SMC,  and  j  is  the  electric  current  density  in  the  core. 

The  copper  loss  is  3.56  (watt)  for  this  PMTG,  and  the  total  resistance  of  coils  is  4.295  (Q).  In  this 
research,  the  skin  effect  is  neglected.  Figure  9  shows  the  profile  of  input  force  in  terms  of  displacement  of 
translator  in  one  reciprocating  cycle,  and  Figure  10  shows  the  load  voltage. 
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Figure  9.  Input  force  versus  displacement 


Figure  10.  Output  voltage 


4.  ANALYTICAL  CALCULATION  OF  LOSSES 

The  losses  are  only  calculated  in  stator  core  and  the  losses  in  the  translator  part  are  negligible.  The 
total  losses  include  the  iron  and  copper  losses,  and  the  iron  loss  is  also  including  the  eddy-current  and 
hysteresis  losses. 
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The  iron  loss  is  calculated  by  using  the  equation  extracted  from  the  catalog  of  SMC  material.  To 
calculate  the  iron  loss,  the  stator  core  is  divided  into  different  parts  with  respect  to  the  flux  density 
distribution.  Since  the  stator  core  is  including  the  yoke  and  tooth  parts,  the  maximum  flux  density  in  each 
part  is  calculated  by  using  the  analytical  method  presented  in  [18].  Figure  11  shows  the  core  partitioning  for 
calculating  the  core  loss  analytically: 


p 

^  FE 


--KjEf  +  KFBl  + 


Blfd^ 


8. 1 X  yc>  X  resistiv  ip  x  1 000 


(9) 


The  comparison  of  results  obtained  through  analytical  method  and  FEM  is  shown  in  Figure  12.  The 
resistance  and  copper  loss  of  coils  can  be  calculated  while  having  the  turn  number,  the  cross  sectional  area, 
and  conductivity  of  coils.  The  obtained  value  for  resistance  is  4.286  (Q),  and  the  copper  loss  is  equal  to  3.6 
(watt).  The  analytically  calculated  loss  including  iron  and  copper  loss  is  equal  to  9.258  (watt). 


Figure  1 1 .  Stator  core  partitioning  for  iron  loss  calculation 


6 

5 

4 


2 

1 

0 


Hysteresis  Eddy  cuiTent  Copper 

loss  loss  loss 


Figure  12.  Loss  comparison  obtained  through  analytical  method  and  FEM 


5.  EQUIVALENT  CIRCUIT  OF  THE  RECIPROCATING  SINGLE-PHASE  PMTG 

Figure  13  shows  the  electric  equivalent  circuit  of  single -phase  reciprocating  PMTG  in  which  Emf 
shows  the  electric  motive  force,  X^is  the  synchronous  reactance,  and  Rg  is  the  stator  resistance: 

+  (10) 

at 


_  .  ..dcp  ..dtp  dz  ^jdcp 

Emf^-N—  =  -N—  X  —  =  -N—  x  v 
dt  dz  dt  dz 


(11) 


where  z  ,  v  and  (p  are  displacement,  velocity  and  magnetic  flux,  respectively.  The  synchronous  reactance  is 
calculated  by  using  the  results  obtained  through  short-circuit  and  open-circuit  tests  [19].  The  obtained  values 
of  synchronous  impedance  and  synchronous  inductance  are  16.87  (Q)  and  112  (mH),  respectively. 
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Figure  13.  Equivalent  circuit  of  single -phase  reciprocating  PMTG 


6.  CONCLUSION 

In  this  paper,  the  studied  reciprocating  PMTG  with  solid  SMC  core  was  simulated  and  analyzed  by 
using  nonlinear  transient  2-D  FEM,  in  which  the  J-A  model  of  magnetic  hysteresis  was  incorporated  and 
then,  the  hysteresis  losses  were  calculated.  Meanwhile,  the  eddy-current  and  copper  losses  were  calculated. 
Due  to  using  the  SMC  core,  the  eddy-current  losses  are  negligible  despite  the  using  of  solid  core.  The  losses 
obtained  through  FEM  agree  well  with  the  corresponding  values  calculated  analytically.  The  electric 
equivalent  circuit  of  the  PMLG  was  also  extracted.  With  respect  to  the  obtained  results  in  this  paper,  the 
using  of  SMC  core  in  free-piston  generators  is  recommended. 
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